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We present an introduction to mass and angular momentum in General Rela- 
tivity. After briefly reviewing energy-momentum for matter fields, first in the 
flat Minkowski case (Special Relativity) and then in curved spacetimes with 
or without symmetries, we focus on the discussion of energy-momentum for 
the gravitational field. We illustrate the difficulties rooted in the Equivalence 
Principle for defining a local energy-momentum density for the gravitational 
field. This leads to the understanding of gravitational energy-momentum and 
angular momentum as non-local observables that make sense, at best, for ex- 
tended domains of spacetime. After introducing Komar quantities associated 
with spacetime symmetries, it is shown how total energy-momentum can be 
unambiguously defined for isolated systems, providing fundamental tests for 
the internal consistency of General Relativity as well as setting the concep- 
tual basis for the understanding of energy loss by gravitational radiation. 
Finally, several attempts to formulate quasi-local notions of mass and an- 
gular momentum associated with extended but finite spacetime domains are 
presented, together with some illustrations of the relations between total and 
quasi-local quantities in the particular context of black hole spacetimes. This 
article is not intended to be a rigorous and exhaustive review of the subject, 
but rather an invitation to the topic for non-experts. In this sense we follow 
essentially the expositions in [1, 2, 3, 4] and refer the reader interested in 
further developments to the existing literature, in particular to the excellent 
and comprehensive review by Szabados [1]. 



1 Issues around the notion of gravitational energy in 
General Relativity 

1.1 Energy- momentum density for matter fields 

Let us first consider mass and angular momentum associated with matter in 
the absence of gravity, in a flat Minkowski spacetime. The density of energy 
and linear momentum associated with a distribution of matter are encoded in 
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the energy- momentum tensor T^j^, corresponding to the Noether current con- 
served under infinitesimal spacctimc translations in a Lagrangian framework. 
This general conservation property, namely d^T^'^ = in inertial Minkowski 
coordinates, plays a key role in our discussion. Indeed, together with the 
presence of symmetries, it permits the introduction of conserved quantities 
or charges. Given a space-like hypersurface S and considering the unit time- 
like vector normal to it, we can define the conserved quantity associated 
with the symmetry and the domain D {<Z S) as 

QD[k^] = I kPT^pU^^ifx , (1) 
Jd 

where ^ (Px denotes the induced volume element in D. The conservation 
of T^n and the characterisation of k^ as a symmetry imply the conservation 
of the vector T^'^k'', i.e. {T^'^k'') = 0. Applying then the Stokes theo- 
rem, it follows the equality between the change in time of Qnlk^] and the 
flux of j^pTP'^ki, through the boundary of D (where j'^^ is the projector on 
D). Minkowski spacetime symmetries are given by Poincare transformations. 
Therefore, we can associate conserved quantities with the infinitesimal gen- 
erators corresponding to translations T^, rotations J^, and boosts (here 
the label a for translation generators runs in {0, 1, 2, 3}, whereas i is a space- 
like index in {1, 2, 3}). In this manner, a 4-momentum Paii*] and an angular 
momentum Ji[D] associated with the distribution of matter in D C S can 
be defined as 

Pa[D]= f T^^Tln^^d?x , Ji[D]= f T^.J'fn''^ d^'x . (2) 
Jd Jd 

More generally, we can combine together the rotation and boost generators 
and into a vector-field- valued antisymmetric matrix Mj^j^j (where = 

^ei'''^M|^.^j and = M|^.j) and write the conserved quantities 

J[ab] [D] = T^.M|;,]n''V7 d^x . (3) 

The mass and (Pauli-Lubanski) spin are constructed as 

m''[D]:=-^-^P,[D]P,[D] , S''[D] := ^\'^>"^''P,[D]J[,4D] , (4) 

in terms of which Poincare Casimirs (invariant under Poincare transforma- 
tions) can be expressed. 

In the non-flat case, (matter) energy-momentum tensor acts as the source 
of gravity through the Einstein equation and, consistently with Bianchi iden- 
tities, satisfies the divergence- free condition analogous to the fiat conservation 
law: 

G^.~%.-^^Rg^, = 87rT^, , V,T-^^ = 0. (5) 
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The same strategy employed in the flat case for deflning physical quantities 
associated with matter, i.e. from conserved cmrents corresponding to some 
symmetry, can be followed in non-flat spacetimes {M,g^u) presenting Killing 
vectors fc**. The vector T^^^k" is conserved, i.e. {T'^^k") = 0, and provides 
a current-density for the conserved quantity Qoik'^] defined by expression (1). 
The physical interpretation of Qoik^] depends of course on the nature of the 
Killing vector k^. Actually, Qnlk^] does not actually depend on the slice S 
in the sense that its value is the same in the domain of dependence of D (this 
precisely corresponds to the conserved nature of this charge). 

In a general spacetime with no symmetries the previous strategy ceases to 
work, and ambiguities in the definition of mass and angular momentum enter 
into scene. One can still calculate the flux of T'^^^'^ for a given vector ^'^ , and 
define the associated quantity Qd[^^]- However, the latter will now depend 
on the slice S and, in addition, its explicit dependence on introduces 
some degree of arbitrariness in the discussion. In this context, given a space- 
like 3-1-1 foliation {St} of the spacetime with time-like normal vector n^, 
the current P^^ := — T'"'n,y can be interpreted as the energy-momentum 
density associated with (Eulerian) observers at rest with respect to Sf That 
is, E := T^'^UfiTii, stands as the matter energy density and := —^^ pT'"^nv 
as the momentum density, where is the induced metric on Et (see Eq. 
(12) below for the complete 3-1-1 decomposition of T^v). In particular, we 
can calculate the matter energy associated with observers n'* over the spatial 
region D by direct integration 



By imposing the dominant energy condition on the matter energy-momentum 
tensor (see section 3.3), the vector —T'^^rii, is future directed and non-space- 
like. Its Lorentzian norm is therefore non-positive and an associated matter 
mass density m can be given as := —P'^Pn = — (— T^''np)(— T'"^no-).gpi/ = 
—p^Pi > 0. The corresponding mass M[D] in the extended region D would 



Note the difference between the construction of M[D] and that of m[D] in 
the Minkowskian case: for the latter one first integrates to obtain the charges 
and then calculates a Minkowskian norm, whereas for constructing M[D] that 
order is reversed; in addition, different metrics are employed in each case (cf. 
section 2.2. in [1]). 

1.2 Problems when defining a gravitational energy-momentum 

In the characterisation of the physical properties of the gravitational field, 
in particular its energy-momentum and angular momentum, we could try to 




(6) 



be 
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follow a similar strategy to that employed for the matter fields. This would 
amount to identify appropriate local densities that would then be integrated 
over finite spacetime regions. However such an approach rapidly meets im- 
portant conceptual difficulties. 

A local (point-like) density of energy associated with the gravitational 
field cannot be defined in General Relativity. Reasons for this can be tracked 
to the Equivalence Principle. Illustrated in a heuristic manner, this principle 
can be used to get rid of the gravitational field on a given point of spacetime. 
Namely, a free falling point-like particle does not feel any gravitational field 
so that, in particular, no gravitational energy density can be identified at 
spacetime points. 

In a Lagrangian setting, these basic conceptual difficulties are reflected 
in the attempts to construct a gravitational energy-momentum tensor, when 
mimicking the methodological steps followed in the matter field case. We can 
write generically the gravitational-matter action as 



S — SeU + 'S'm — „ 

(8) 

where Seh denotes the Einstein-Hilbert action and in the matter La- 
grangian Lm account for the matter fields. The symmetric energy-momentum 
for matter is obtained from the variation of the matter action with respect 
to the metric 

-2_6S^ 

whereas the field equations for the matter fields follow from the variation 
with respect to the matter fields On the contrary, the gravitational ac- 
tion ^EH only depends on the gravitational field, since any further background 
structure would be precluded by diffeomorphism invariance (a feature closely 
tied to the physical Equivalence Principle). Einstein equation for the grav- 
itational field follows from the variation of the total action with respect to 
the metric field g^j^, with no gravitational analogue of the symmetric matter 
energy-momentum tensor Tf^^. Attempts to construct a symmetric energy- 
momentum tensor for the gravitational field either recover the Einstein tensor 
or can only be related to higher-order gravitational energy-momentum 
objects, such as the Bel- Robinson tensor (see e.g. [5]). Again, the absence of a 
tensorial (i.e. point-like geometric) quantity representing energy-momentum 
for the gravitational field is consistent with, and actually a consequence of, 
the Equivalence Principle. 

The natural interpretation of the symmetric matter energy-momentum 
tensor T^,^ as introduccid in Eq. (9) is that of the current- source for the 
gravitational field, obtained as a conserved current associated with spacetime 
translations. Alternative, in terms of the Noether theorem [6] it is natural to 
introduce a (non-symmetric) canonical energy-momentum tensor for matter 
from which a symmetric one can be constructed through the Belinfante- 
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Rosenfeld procedure [7, 8, 9]. The application of this construction to the 
gravitational field naturally leads to the disciission of gravitational energy- 
momentum pseudo-tensors [1]. The underlying idea consists in decomposing 
the Einstein tensor G^i, into a part that can be identified with the energy- 
momentum and a second piece that can be expressed in terms of a pseudo- 
potential. That is [10] 

G/:=-87rt/ + ^^5^(if/^) , (10) 

where i^i^ is the gravitational energy- momentum pseudo-tensor and .ff^i,^ is 
the superpotential. Einstein equation is then written as 

aA(^f/^) = 167^^/=^(^/ + V)=:167rr/ . (11) 

Objects tfj,'^ and H^"^^ are not tensorial quantities. This means that their value 
at a given spacetime point is not a well-defined notion. Moreover, their very 
definition needs the introduction of some additional background structure and 
some choice of preferred coordinates is naturally involved. Different pseudo- 
tensors exist in the literature, e.g. those introduced by Einstein, Papapetrou, 
Bergmann, Landau and Lifshitz, Moller or Weinberg (e.g. see references in 
[10]). 

As an alternative to the pseudo-tensor approach, there also exist attempts 
in the literature aiming at constructing truly tensorial energy-momentum 
quantities. However they also involve the introduction of some additional 
structure, either in the form of a background object or by fixing a gauge in 
some given formulation of General Relativity (cf. comments on the tetrad 
formalism approach in [1]). 

Non-local character of gravitational energy 

As illustrated above, crucial conceptual and practical caveats are involved in 
the association of energy and angular momentum with the gravitational field. 
For these reasons, one might legitimately consider gravitational energy and 
angular momentum in General Relativity as intrinsically meaningless notions 
in generic situations, in such a way that the effort to derive explicit general 
local expressions actually represents an ill-defined problem (cf. remarks in [11] 
referring to the quest for a local expression of energy in General Relativity) . 
Having said this and after accepting the non-existence of a local (point-like) 
notion of energy density for the gravitational field, one may also consider 
gravitational energy-momentum and angular momentum as notions intrinsi- 
cally associated with extended domains of the spacetime and then look for 
restricted settings or appropriate limits where they can be properly defined. 

In fact, making a sense of the energy and angular momentum for the 
gravitational field in given regions of spacetime is extremely important in 
different contexts of gravitational physics, as it can be illustrated with exam- 
ples coming from mathematical relativity, black hole physics, lines of research 
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in Quantum Gravity, or relativistic astrophysics. Prom a structural point of 
view, having a well-defined mass positivity result is crucial for the internal 
consistency of the theory, as well as for the discussion of the solutions sta- 
bility. Moreover, the possibility of introducing appropriate positive-definite 
(energy) quantities is often a key step in different developments in mathemat- 
ical relativity, in particular when using variational principles. In the study of 
the physical picture of black holes, appropriate notions of mass and angular 
momentum are employed. In particular, they play a key role in the formu- 
lation of black hole thermodynamics (e.g. [12]), a cornerstone in different 
approaches to Quantum Gravity. In the context of relativistic astrophysics 
and numerical relativity, the study of relativistic binary mergers, gravita- 
tional collapse and the associated generation/propagation of gravitational 
radiation also requires appropriate notions of energy and angular momentum 
(sec e.g. [13] for a further discussion on the intersection between numerical 
and mathematical relativity). 

Once the non-local nature of the gravitational energy-momentum and an- 
gular momentum is realised, the conceptual challenge is translated into the 
manner of determining the appropriate physical parameters associated with 
the gravitational field in an extended region of spacetime. An unambiguous 
answer has been given in the case of the total mass of an isolated system. 
However, the situation is much less clear in the case of extended but finite 
spacetime domains. In a broad sense, existing attempts either enforce some 
additional structure that restricts the study to an appropriate subset of the 
solution space of General Relativity, or alternatively they look for a genuinely 
geometric characterisation aiming at fulfilling some expected physical require- 
ments. In this article wc present an overview of some of the relevant existing 
attempts and illustrate the kind of additional structures they involve. 

1.3 Notation 

Before proceeding further, we set the notation, some of whose elements have 
already been anticipated above. The signature of spacetime {M.,g^i,) is cho- 
sen to be diag[ — 1,1,1,1] and Greek letters are used for spacetime indices 
in {0, 1, 2, 3}. We denote the Levi-Civita connection by V;^ and the volume 
element by = dx^ A dx^ A dx'^ A dx^. We make G = c = 1 throughout. 

3+1 decompositions 

In our presentation of the subject, 3-1-1 foliations of spacetime {M.^g^^) by 
space- like 3-slices {St} will play an important role. Given a height- function 
t, the time- like unit normal to St will be denoted by n'' and the 3+1 de- 
composition of the evolution vector field by = Nn^ + where N is the 
lapse function and is the shift vector. The induced metric on the space- 
like 3-slice St is expressed as 7^1/ = g^v + ni^n„, with Z)^ the associated 
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Levi-Civita connection and volume element = yfydx^ A dx^ A dx^, so that 
^e^^p = "-"^^ecr/xi/p- The extrinsic curvature of (I7t,7^y) in {A4,g^^) is defined 
as Kij,i, := —\Cn''ifj,v = — 7/i''Vpny. The 3+1 decomposition of the (matter) 
stress-energy tensor, in terms of an Eulerian observer n'' in rest with respect 
to the foliation {St}-, is 

Tu,i, = E n^n„ + p^^n^) + S^v , (12) 

where the matter energy and momentum densities are given by E := T/^vn^n'^ 
andp^ := —T,yf/n""j''^^, respectively, whereas the matter stress tensor is S^^v ■= 
Tpa'y^ ^'y'^ Latin indices running in {1, 2, 3} will be employed in expressions 
only involving objects intrinsic to space-like St slices. 

Closed 2-surfaces 

Closed 2-surfaces S, namely topological spheres in our discussion, will also be 
relevant in the following. The normal bundle T-^S can be spanned by a time- 
like unit vector field and a space- like unit vector field s^, that we choose to 
satisfy the orthogonality condition n^s^ = 0. When considering S as embed- 
ded in a space-like 3-surface S, n'^ can be identified with the time-like normal 
to S and with the normal to S tangent to S. In the generic case, and 

can be defined up to a boost transformation: n'^ = cosh{r])n^ + sinh(ry)s'* 
and s'^ = sinh(77)n^ + cosh(77)s'^, with rj a real parameter. Alternatively, one 
can span Tj^S at p G 5 in terms of the null normals defined by the inter- 
section between the normal plane to S and the light-cone at the spacctime 
point p. The directions defined by the outgoing i'^ and the ingoing null 
normals (satisfying k^t^ = —1) are uniquely determined, though it remains a 
boost- normalization freedom: = k'^ = j'k^- The induced metric on 

S is given by: q^^ = g^^ + k^t^+l^k^ = gtiv + n^,n^ - s^s„ = 7^1, - s^s^, the 
latter expression applying when S is embedded in (£",7^^). The Levi-Civita 
connection associated with q^y will be denoted by ^D^ and the volume ele- 
ment by \ = y^dx^ A dx^, i.e. ^e^^ — s'^'^f-paftu- When integrating tensors 
on S with components normal to the sphere, it is convenient to express the 
volume element as dSp,^ = {sp-n^ — npSy)^/qd^x (this is just a convenient 
manner of re-expressing \pup(T for integrating over <S after a contraction with 
the appropriate tensor; cf. for example Eq.(13)). 

The second fundamental tensor oi (<S, qij,^) in {M,gfj,u) is defined as /C^^, := 

q^ff ,y pl^a^ that can be expressed as = n'^O^pJ + s"0^!] = k'^O^^} + 
i'^O'^p} , where the deformMion tensor O^^J associated with a vector normal 

(v) 

to <S is defined as OpJ = q'' ^q'^ pVc,. We set a specific notation for the cases 

corresponding to and n**, namely ff^y := 0^*^, the extrinsic curvature of 

{S,q^u) inside a 3-slice {S,')^^), and := -0^"^ 

Information about the extrinsic curvature of {S^qf^v) in (Ai^gp,!,) is com- 
pleted by the normal fundamental forms associated with normal vectors v'^ . In 



8 



Jose Luis Jaramillo and Eric Gourgoulhon 



particular, we define the 1-form Oj;^^ := k'^q'^ „lp. This form is not invari- 
ant under a boost transformation, and transforms as + '^Df^lnf in 

the notation above. Other normal fundamental forms can be defined in terms 
of normals k'^,n'^ and s'', but they are all related up to total derivatives. 

2 Spacetimes with Killing vectors: Komar quantities 

As commented above, some additional structure is needed to introduce mean- 
ingful notions of gravitational energy and angular momentum. Let us first 
consider spacetimes admitting isometries. This represents the most straight- 
forward generalization of the definition of physical parameters as conserved 
quantities under existing symmetries. Requiring the presence of Killing vec- 
tors represents our first example of the enforcement of an additional structure 
on the considered spacetime. 

Given a Killing vector field /c^ in the spacetime {M.,g^u) and S a space- 
like closed 2-surface, let us define the Komar quantity [14] A;k as 

k^~-^j^V^k''dS^, , (13) 

(see previous section for the notation dS^v for the volume element on S). 
Let us consider <S as embedded in a space-like 3-slice E and let us take a 
second closed 2-surface S' such that either S' is completely contained in <S or 
vice- versa, and let us denote by V the region in S contained between S and 
S'. The previously defined Komar quantity kx is then conserved in the sense 
that its value does not depend on the chosen 2-surface as long as no matter 
is present in the intermediate region V 

fcg = (t^, - ^Tg^,^ n'^k-'^d^x + fcg' , (14) 

where T = T^^g'"'. 

Remark 1. Two important points must be stressed: a) the definition of fcx 
is geometric and therefore coordinate independent, and b) A;k is associated 
with a closed 2-surface with no need to refer to any particular embedding in 
a 3-slice S (in the discussion above the latter has been only introduced for 
pedagogical reasons) . 

2.1 Komar mass 

Stationary spacetimes admit a time-like Killing vector field fc**. The associated 
conserved Komar quantity is known as the Komar mass 



Mk :=-^£v^rd5^. 



(15) 
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This represents our first notion of mass in General Relativity. It is instructive 

to write the Komar mass in terms of 3+1 quantities. Given a 3-shcing {St} 
and choosing the evolution vector = Nn^ + to coincide with the time- 
hke Killing symmetry, we find 

M^ = ^i {s'DiN - Kijs'/3^) ^d^x . (16) 



2.2 Komar angular momentum 

Let us consider now an axisymmetric spacetime, where the axial Killing vector 
is denoted by That is, is a space- like Killing vector whose action on M 
has compact orbits, two stationary points (the poles), and is normalized so 
that its natural afhne parameter takes values in [0,27r). The Komar angular 
momentum is defined as 

Jk:=y^£ V^rdS^^ . (17) 

Note (apart from the sign choice) the factor 1/2 with respect to the Komar 
quantity ^k, known as the Komar anomalous factor (it can be explained in 
the context of a bimetric formalism by writing the conserved quantities in 
terms of an Einstein energy-momentum flux density that can be expressed as 
the sum of half the Komar contribution plus a second term: in the angular 
momentum case this second piece vanishes, whereas for the mass case it equals 
half the Komar term; cf. [15]). Adopting a 3-slicing adapted to axisymmetry, 
i.e. n^<pi_i = 0, we have: 

J^ = ^j^Kijs'(l>'^dPx=j-j^f2\^^(^^^Sx . (18) 



3 Total mass of Isolated Systems in General Relativity 
3.1 Asymptotic Flatness characterisation of Isolated Systems 

The characterisation of an isolated system in General Relativity aims at cap- 
turing the idea that spacetime becomes flat when we move sufficiently far 
from the system, so that spacetime approaches that of Minkowski. However, 
the very notion of far away becomes problematic due to the absence of an a 
priori background spacetime. In addition, we must consider different kinds of 
infinities, since we can move away from the system in space-like and also in 
null directions. Different strategies exist in the literature for the formaliza- 
tion of this asymptotic flatness idea, and not all of them are mathematically 
equivalent. Traditional approaches attempt to specify the adequate fall-off 
conditions of the curvature in appropriate coordinate systems at infinity. 
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These approaches have the advantage of embodying the weakest versions 
of asymptotic flatness. Wc will illustrate their use in the discussion of spa- 
tial infinity in section 3.2. Ifowever, the use of coordinate expressions in 
this strategy also introduces the need of verifying the intrinsic nature of the 
obtained results, something that it is not always straightforward. For this 
reason, a geometric manner of describing asymptotic flatness is also desir- 
able, without relying on specific coordinates. This has led to the conformal 
compactification picture, where infinity is brought to a finite distance by an 
appropriate spacetime conformal transformation. More concretely, one works 
with an unphysical spacetime {M,gij,v) with boundary, such that the physical 
spacetime {M,g^„) is conformally equivalent to the interior of {M,g^^), i.e. 
g^i, = fi^Qfjtv Infinity is captured by the boundary dM. and is characterised 
by the vanishing of the conformal factor, i? = 0. The whole picture is inspired 
in the structure of the conformal compactification of Minkowski spacetime. 
The conformal boundary is the union of different pieces, which are classi- 
fied according to the metric-type of the geodesies reaching their points. This 
defines (past and future) null infinity spatial infinity i° and (past and 
future) time-like infinity i^, i.e. dM. = J^^Ui^Ui^ . The conformal spacetime 
is represented in the so-called Carter-Penrose diagram. Fall-off conditions for 
the characterisation of asymptotic fiatness arc substituted by differentiability 
conditions on the fields at null and spatial infinity (isolated systems do not 
require flatness conditions on time- like infinity) . Null infinity was introduced 
in the conformal picture by Penrose [16, 17], the discussion of asymptotic 
fiatness at spatial infinity was developed by Geroch [18] and a unified treat- 
ment was presented in [19, 20] (see also [21, 22]). We will briefiy illustrate the 
different approaches to asymptotic flatness in the following sections, but we 
refer the reader to the existing bibliography (e.g. [23, 4]) for further details. 

3.2 Asymptotic Euclidean slices 

The following two sections are devoted to the discussion of conserved quanti- 
ties at spatial infinity, but they also illustrate the coordinate-based approach 
to asymptotic flatness. A slice S endowed with a space-like 3-metric ^ij is 
asymptotically Euclidean (flat), if there exists a Riemannian background met- 
ric fij such that: 

i) fij is flat, except possibly on a compact domain D of E. 

ii) There exists a coordinate system (x*) = {x,y,z) such that outside D, 
fij = diag(l,l,l) {Cartesian-type coordinates) and the variable r := 
y^a?2+j/2 + z2 (,.^]^ ^g^jjg arbitrarily large values on E. 

iii) When r — > -l-oo 

Kij = 0{r-^) 



dx^ 



0{r-') 
= 0(r-3) 



(19) 
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Given an asymptotically flat spacetime foliated by asymptotically Euclidean 
slices {-S'tji spatial infinity is deflned by r ^ +00 and denoted as i^. 

Asymptotic symmetries at spatial infinity 

As commented in the discussion of the Komar quantities, the existence of 
symmetries provides a natural manner of defining physical parameters as 
conserved quantities. In the context of spatial infinity, the spacetime diffeo- 
morphisms preserving the asymptotic Euclidean structure (19) are referred 
to as asymptotic symmetries. Asymptotic symmetries close a Lie group. 
Since the spacetime is asymptotically flat, one would expect this group to 
be isomorphic to the Poincare group. However, the set of diffeomorphisms 
(x^) = {t,x^) —7- {x'^) — {t',x'^) preserving conditions (19) is given by 

a;"^ = yl^x" + c^(0,(^) + O(r-i) , (20) 

where A^^^ is a Lorentz matrix and the c'^'s are four functions of the an- 
gles {d,(p) related to coordinates (x*) = {x,y,z) by the standard spherical 
formulae: x — r sind cosip,y = r sinO simp, z = rcos9. This group indeed 
contains the Poincare symmetry, but it is actually much larger due to the 
presence of angle- dependent translations. The latter arc known as supertrans- 
lations and arc defined by c''(0, ^ const and yl^,^ = 5^^, in the group rep- 
resentation (20). The corresponding abstract infinite-dimensional symmetry 
preserving the structure of spatial infinity {Spi) is referred to as the Spi group 
[19, 20]. The existence of this (infinite-dimensional) Lie structure of asymp- 
totic symmetries has implications in the definition of a global physical mass, 
linear and angular momentum at spatial infinity (see below). 

3.3 ADM quantities 

Hamiltonian techniques are particularly powerful for the systematic study 

of physical parameters, considered as conserved quantities under symmetries 
acting as canonical transformations in the solution (phase) space of a theory. 
In this sense, the Hamiltonian formulation of General Relativity provides a 
natural framework for the discussion of global quantities at spatial infinity. 
This was the original approach adopted by Arnowitt, Deser and Misner in 
[24] and we outline here the basic steps. 

First, a variational problem for the class of spacetimes we are considering 
must be set. For a correct formulation we need to specify: a) the dynamical 
fields we are varying, b) the domain V over which these fields are varied 
together with the prescribed value of their variations at the boundary 9V, and 
c) the action functional S compatible with the field equations. As integration 
domain V we consider the region bounded by two space-like 3-slices St^ 
and St2 s-nd an outer time-like tube B. Sti and can be seen as part 
of a 3-slicing {St} with metric and extrinsic curvature given by {jij,K''-^), 
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whereas B has {xi^v^P^'^) as induced metric and extrinsic curvature. That 
iSj X^^l' — Qfiiy ^ u^Ui, and P^^, ~ ^Ip!''^ pUy^ where is the unit space- hke 
normal to B. The dynamical field whose variation we consider is the spacetime 
metric 5^1,, under boundary conditions ^gp.v\dv = (note that we impose 
nothing on variations of the derivatives of g^,y). The appropriate gravitational 
Einstein- Hilbert action then reads (cf. for example [3]; the discussion has a 
straightforward extension to incorporate matter) 



(21) 



^ = ll //^^^'^ + ^ { - - ""^^^''^ 

+ [ {K-Ko)^d'x+ [ {P-Po)V^d^x\ , 
Jst^ Jb J 

where K and P are the traces of the extrinsic curvatures of the hypersurfaces 
Sti and B, respectively, as embedded in {M,gpu)- The subindex corre- 
sponds to their extrinsic curvatures as embedded in [M^rj^y). The boundary 
term guarantees the well-posedness of the variational principle, i.e. the func- 
tional differentiability of the action and the recovery of the correct Einstein 
field equation, under the assumed boundary conditions for the dynamical 
fields. 

Making use of the 3+1 fields decompositions, and considering the in- 
tersections St := S n St between space-like 3-slices Et and the time-like 
hypersurface B, we can express the action (21) as 

r {I N{^R + KijK'i - K^) y/^d'x + 2 I {H-Ho)^ d^xl dt 

(22) 

where H and Hq denote the trace of the extrinsic curvature of the 2-surface St 
as embedded in {St,jij) and {St,fij), respectively. The Lagrangian density 
L can be read from the form of the action (22). The 3-metric jij plays the 
role of the dynamical variable and the dependence of L on 7^ follows from 
the explicit expression of the extrinsic curvature Kij in terms of the lapse 
and the shift, that is 

= ^(7ife-Di/3'+7,fcA^'=-7ii) • (23) 

In particular no derivatives of N and appear in (22), indicating that the 
lapse function and the shift vector are not dynamical variables. The Hamil- 

tonian description is obtained by performing a Legendre transformation from 
variables (7ij,7y) to canonical ones (7^,77*-'), where 



The Hamiltonian density "H is then given by 
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'H=T,ij^..-L , (25) 

and the Hamiltonian follows from an integration over a 3-slice, resulting in 
(cf. [2, 3] for details) 



(26) 



-2£ [N{H ~ Ho) - l3'{K,j ~ Kj,j)sq ^(fx^ , 

where 

Ci:=DjK\-DiK . (27) 

Functionals Cq and Cj vanish on solutions of the Einstein equation (in vac- 
uum). More specifically, equations Cq = and Ci = respectively represent 
the Hamiltonian and momentum constraints of General Relativity, corre- 
sponding to the contraction of the Einstein equation (5) with n''. Prom a 
geometric point of view, they arc referred to as the Gauss-Codazzi relations 
and represent conditions for the embedding of (i3't,7y) as a submanifold of 
a spacetime {M.,gi^,^) with vanishing n^'G^j/. The evaluation of the gravita- 
tional Hamiltonian (26) on solutions to the Einstein equation yields 



/^solution = £ [^(^ - Ho) - P\Ki, - Kji,)sq ^Cfx 



(28) 



Remark 2. Note that in the absence of boundaries the gravitational Hamil- 
tonian vanishes on physical solutions. This is a feature of diffeomorphism 
invariant theories [25] and reflects the fact that the Hamiltonian, considered 
as the generator of a canonical transformation, docs not move points in the 
solution space of the theory. In other words, it is a generator of gauge trans- 
formations, something consistent with the interpretation of the Hamiltonian 
as the generator of diffcomorphisms. Note also that the situation changes 
in the presence of boundaries, where diffcomorphisms not preserving bound- 
ary conditions do not correspond to gauge transformations, indicating the 
presence of residual degrees of freedom (this is of relevance, for instance, in 
certain aspects of the quantum theory). 



ADM energy 

We focus on solutions corresponding to isolated systems and consider 3-slices 
St that are asymptotically Euclidean in the sense of conditions (19) (we refer 
the reader to [26] for a discussion of the total energy in cosmological asymp- 
totically Anti-de Sitter spacetimes and to [27] for its discussion in higher 
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curvature gravity theories). We choose the lapse and the shift so that the 

evolution vector is associated with some asymptotically incrtial observer 
for which = 1 and = at spatial infinity. In particular, this flow vec- 
tor generates asymptotic time translations that, in this asymptotically flat 
context, constitute actual (asymptotic) symmetries. Conserved quantities un- 
der time translations have the physical meaning of an energy. In the present 
case, the conserved quantity is referred to as the ADM energy. The latter is 
obtained from expression (28) by making = 1 and /3* = and taking the 
limit to spatial inflnity, namely r — >■ oo in the well-defined sense of section 
3.2. That is 

-Eadm---^ lim <j>{H~Ho)^d^x. (29) 

»7r 5(t,r^oo) J St 

This ADM energy represents the total energy contained in the slice Using 
the explicit expression of the extrinsic curvature in terms of metric compo- 
nents, the ADM energy can be written as 

EAJiM = ^ lim / [V^lij-Vi{f''^^kl)\8'y/q(fx , (30) 

where Vi stands for the connection associated with the metric and, con- 
sistently with notation in section 1.3, s' corresponds to the unit normal to 
St tangent to St and oriented towards the exterior of St (note that when 
r — > oo the normalization with respect to 7^ and fij are equivalent). In par- 
ticular, if we use the Cartesian-like coordinates employed in (19) we recover 
the standard form (see e.g, [4]) 

Remark 3. We note that asymptotic flatness conditions (19) guarantee the 
finite value of the integral since the O(r^) part of the measure y/qcPx is com- 
pensated by the 0(r^^) parts of d^ij/dx^ and d^jj/dx'^. It is very important 
to point out that finiteness of the ADM energy relies on the subtraction of 
the reference value Hq in Eq. (29). 

Conformal decomposition expression of the ADM energy. 

A useful expression for the ADM energy in certain formulations of the Ein- 
stein equation is given in terms of a conformal decomposition of the 3-metric 

li3='^%i ■ (32) 

Choosing the representative 7^ of the conformal class by the unimodular 
condition det(7jj) = det(/ij) = 1, conditions (19) translate into 
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liJ=h+Oir-') , |^ = 0(r-^) , (33) 
for the conformal factor and the conformal metric. Then it follows [2] 

Eadm = -^ lim I s'{vi<P-\v3%A^(fx . (34) 

in S(t,r^oo) J St \ ^ J 

Note that, whereas in the time-symmetric {K^^, = 0) conformally flat case 
the Komar mass is given in terms of the monopolar term in the asymptotic 
expansion of the (adapted) lapse, the ADM energy is given by the monopolar 
term in ip (the latter holds more generally under a vanishing Dirac-like gauge 

condition on V^^/ij). 

Example 1 (Newtonian limit). As an application of expression (34) we check 
that the ADM energy recovers the standard result in the Newtonian limit. 
For this we assume that the gravitational field is weak and static. In this 
setting it is always possible to find a coordinate system [x'^) = [x^ = ct,x^) 
such that the metric components take the form 

-rfr^ = g^^dx^'dx" = - (1 + 2#) dt^ + (1 - 2#) fij dx'dx^ , (35) 

where again is the flat Euclidean metric in the 3-dimensional slice and 
^ is the Newtonian gravitational potential, solution of the Poisson equation 
A<1> = A-Kp where p is the mass density (we recall that we use units in which 
the Newton's gravitational constant G and the light velocity c are unity). 
Then, using ^ = {I - 2^)^* w 1 - i^, Eq. (34) translates into 



Ekt,m = ^ lim / s'Vi^y^d^x = ^ [ A^./fd^x . 



(36) 



where in the second step we have assumed that Et has the topology of M.'^ 
and have applied the Gauss-Ostrogradsky theorem (with A = ViV^). Using 
now that is a solution of the Poisson equation, we can write 

^ADM= / p^fd?x , (37) 



and we recover the standard expression for the total mass of the system at 
the Newtonian limit (as it will be seen in next section, in a non-boosted slice 
like this, mass is directly given by the energy expression). 

ADM 4-momentum. ADM mass 

ADM linear momentum 



Linear momentum corresponds to the conserved quantity associated with 
an invariance under spatial translations. In the asymptotically flat case, the 
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ADM momentum is associated with space translations preserving the fall-off 
conditions (19) expressed in terms of the Cartesian-type coordinates (a;*). 
Given one of such coordinate systems, the three vectors (9j)jg{i^2,3} repre- 
sent asymptotic symmetries generating asymptotic spatial translations that 
correspond to a choice N = and = <5] in the evolution vector t^. 

Substituting these values for the lapse and shift in the Hamiltonian expres- 
sion evaluated on solutions (28), we obtain the conserved quantity under the 
infinitesimal translation 9,: 

Pi:=^ lim <f {K,k-K^^k) s'^Vqcfx . (38) 

Remark 4- Asymptotic fall-off conditions (19) guarantee the finiteness of ex- 
pression (38) for Pi. 

The ADM momentum associated with the hypersurface Ut is defined as the 
linear form (Pj) = {Pi,P2,P3). Its components actually transform as those 
of a linear form under changes of Cartesian coordinates (a;*) — > (a;'*) which 
asymptotically correspond to a rotation and/or a translation. For discussing 
transformations under the full Poincare group, we must introduce the ADM 
4-momentum defined as 

{P^''^):={-Eat,m,PuP2,P3) . (39) 

Under a coordinate change (a;'') = (i, a;*) {x'^) = {t',x'^) which preserves 

the asymptotic conditions (19), i.e. any coordinate change of the form (20), 
components p^^M transform under the vector linear representation of the 
Lorentz group 

as first shown by Arnowitt, Deser and Misner in [24]. Therefore (P^^^) can 
be seen as a linear form acting on vectors at spatial infinity and is called 
the ADM ^-momentum. 

ADM mass 

Having introduced the ADM 4-momentum, its Minkowskian length provides 
a notion of mass. The ADM mass is therefore defined as: 

MIdm := -^;f°''^^DM , ^ADM = V^^iDM-^i^' • (41) 

Remark 5. In the literature, references are found where the term ADM mass 
actually refers to this length of the ADM 4-momentum and other references 
where it refers to its time component, that we have named here as the ADM 
energy. These differences somehow reflect traditional usages in Special Rel- 
ativity where the term mass is sometimes reserved to refer to the Poincare 
invariant (rest-mass) quantity, and in other occasions is used to denote the 
boost-dependent time component of the energy-momentum. 
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The ADM mass is a time independent quantity. Time evolution is generated 

by the Hamiltonian in expression (26). The time variation of a given quantity 
F defined on the phase space is expressed as the sum of its Poisson bracket 
with the Hamiltonian (accounting for the implicit time dependence through 
the time variation of the phase space variables) and the partial derivative 
of F with respect to time. Since in expression (26) there is no explicit time 
dependence, constancy of the ADM mass follows: 



As a consequence of this, the ADM mass is a property of the whole (asymp- 
totically flat) spacetime. 

Remark 6 (Relation between ADM and Komar masses). Komar mass is de- 
fined only in the presence of a time-like Killing vector fc^ (more generally, cf. 
[28] for an early critical account of its physical properties). However, in the 
asymptotically flat case wc can discuss the relation between the ADM energy 
and the Komar mass associated with an asymptotic inertial observer. Though 
the relation is not straightforward from explicit expressions (18) and (30), it 
can be shown [29, 21] that, for any foliation {St} such that the associated 
unit normal coincides with the time-like Killing vector k^^ at infinity (i.e. 
A'' ^ 1 and 0) we have 



As a practical application, this relation has been used as a quasi-equilibrium 
condition in the construction of initial data for compact objects in quasi- 
circular orbits (e.g. [30]). 

Positivity of the ADM mass. 

One of the most important results in General Relativity is the proof of the 

positivity of the ADM mass under appropriate energy conditions for the mat- 
ter energy-momentum tensor. This is important first on conceptual grounds, 
since it represents a crucial test of the internal consistency of the theory. A vi- 
olation of this result would evidence an essential instability of the solutions of 
the theory. It is also relevant on a practical level, since this theorem (and/or 
related results) pervades the everyday practice of (mathematical) relativists. 

The theorem states that, under the dominant energy condition, the ADM 
mass cannot be negative, i.e. Madm > 0. Moreover, Madm = if and only if 
the spacetime is Minkowski. This result was first obtained by Schoen and Yau 
[31, 32] and then recovered using spinorial techniques by Witten [33] (sec in 
this sense [34] for a previous related mass positivity result in supergravity). 

The dominant energy condition essentially states that the local energy 
measured by a causal observer is always positive, and that the flow of energy 
associated with this observer cannot travel faster than light. More precisely. 



—Madm = . 



(42) 



Mk = Madm 



(43) 
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given a future-directed time-like vector u**, this conditions states that the 
vector —T^yV^ is a futurc-oricntcd causal vector. Vector —T^^.v'^ represents 
the energy-momentum 4-current density as seen by the observer associated 
with v^, in an analogous decomposition to that in (12). From the dominant 
energy condition it follows E Tf^i.v^^v'^ > 0, i.e. the local density cannot 
be negative {weak energy condition) and, more generally, E > \/P^Pi. 

ADM angular momentum 

Pushing forward the strategy followed for defining the ADM mass and linear 
momentum, one would attempt to introduce total angular momentum as the 
conserved qiiantity associated with rotations at spatial infinity. More specif- 
ically, in the Cartesian-type coordinates used for characterising asymptoti- 
cally Euclidean slices (19), infinitesimal generators (<;^i)je{i,2,3} for rotations 
around the three spatial axes are 

0x = -zdy + ydz , (l)y = -xd^ + zd^ , = -yd^ + xdy , (44) 

which constitute Killing symmetries of the asymptotically flat metric. When 

using the associated lapse functions and shift vectors in the Hamiltonian ex- 
pression (28), namely N = and = the following three quantities 
result 



However, the interpretation of J, as the components of an angular momentum 
faces two problems: 

1. First, asymptotic fall-off conditions (19) are not sufficient to guarantee 
the finiteness of expressions (45). 

2. Second, in contrast with the linear momentum case, the quantity (Jj) = 
(Ji, J2, J3) does not transform appropriately under transformations (20) 
preserving (19). This can be tracked to the existence of supertranslations. 
In particular, the so-deflned angular-momentum vector ( Jj) depends non- 
covariantly on the particular coordinates we have chosen. 

For this reason, it is not appropriate to refer to an ADM angular momentum 
in the same sense that we use the ADM term for mass and linear momen- 
tum quantities. A manner of removing the above-commented ambiguities 
consists in identifying an appropriate subclass of Cartesian-type coordinates 
where, first, the Jj components are finite and, second, they transform as the 
components of a linear form. Among the different strategies proposed in the 
literature, we comment here on the one proposed by York [35] in terms of 
further conditions on the conformal metric jij introduced in (32) and the 
trace of the extrinsic curvature K. Namely 




iG {1,2,3} . (45) 
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(46) 



representing asymptotic gauge conditions. That is, they actually impose re- 
strictions on the choice of coordinates but not on the geometric properties of 
spacetime at spatial infinity. First condition in (46) is known as the quasi- 
isotropic gauge, whereas the second one is referred to as the asymptotic max- 
imal gauge. 

Remark 7. Note that, in contrast with the total angular momentum defined 
at spatial infinity, no ambiguity shows up in the definition of the Komar 
angular momentum in Eq. (17). 

3.4 Bondi energy and linear momentum 

We could introduce Bondi (or Trautman-Bondi-Sachs) energy at null infinity 
following the same approach we have employed for the ADM energy, i.e. by 
taking the appropriate limit of (28) with N = 1 and = 0. In the present 
case, instead of keeping t constant and making r — oo as we did in (29), 
we should introduce retarded and advanced time coordinates (respectively, 
u = t — r and v = t + r) and consider the limit 



The full discussion of this limit would require the introduction of the ap- 
propriate fall-off conditions for the metric components in a special class of 
coordinate system adapted to null infinity (Bondi coordinates). This is in the 

spirit of the original discussion on the energy flux of gravitational radiation 
from an isolated system by Bondi, Van der Burg and Mctzncr [36], and Sachs 
[37]. However, aiming at providing some fiavour of the geometric approach 
to asymptotic flatness, we rather outline here a discussion in the setting of 
the conformal compactification approach. 

Null infinity 

A smooth spacetime (A^, (?) is asymptotically simple [38] (see e.g. also [23]) if 
there exists another (unphysical) smooth Lorentz manifold {A4,g) such that: 

i) Ai is an open submanifold of M with (smooth) boundary dAi. 

ii) There is a smooth scalar field H on M, such that: il > 0, g^v = ^^giiv 



An asymptotically simple spacetime is asymptotically flat (at null infinity) if, 

in addition, Einstein vacuum equation is satisfied in a neighbourhood of dM. 
(or the energy-momentum decreases sufficiently fast in the matter case). In 




(47) 



on M, and 12 = 0, dfj^fi ^ on dM. 
iii) Every null geodesic in M. begins and ends on dM.. 
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this case the boundary dM. consists, at least, of a null hypersurface with two 
connected components ^ = U each one with topology 5^ x M (note 
that in Minkowski dM also contains the points Boundaries and 

represent past and future null infinity, respectively. 

Symmetries at null infinity 

In order to characterise a vector ^'^ in A4 as an infinitesimal asymptotic 
symmetry at (future) null infinity wc must assess the vanishing of £^(7^,^ 
as one gets to For this, wc require first that considered as a vector 
field in the unphysical spacetime (i.e. under the immersion of A4 into A4), 
can be smoothly extended to Then is characterised as an asymptotic 
symmetry by demanding that Q^L^g^i, can also be smoothly extended to 
J^^ and vanishes there, that is 

Two vector fields ^'^ and ^''^ are considered to generate the same infinites- 
imal asymptotic symmetry if their extensions to coincide. The equiva- 
lence class of such vector fields, that we will still denote by generates the 
asymptotic symmetry group at This is known as the Bondi-Metzner- 
Sachs (BMS) group and is universal in the sense that it is same for every 
asymptotically flat spacetime. The BMS group is infinite-dimensional, as it 
was the case of the Spi group at spatial infinity. It does not only contain the 
Poincare group, but actually is a semi-direct product of the Lorentz group 
and the infinite-dimensional group of angle dependent supertranslations (see 
details in e.g. [4]). The key point for the present discussion is that it possesses 
a unique canonical set of asymptotic 4-translations characterised as the only 
4-parameter subgroup of the supertranslations that is a normal subgroup of 
the BMS group. This leads us to the Bondi-Sachs 4-momentum. 

Bondi-Sachs 4-momentum 

As mentioned above, the original introduction of the Bondi energy was based 
in the identification of certain expansion coefficients in the line element of ra- 
diative spacetimes in adapted (Bondi) coordinates [36]. A Hamiltonian anal- 
ysis, counterpart of the approach adopted in section 3.3 for introducing the 
ADM mass, can be found in [39]. Here we rather follow a construction based 
on the Komar mass expression. Though Eq. (13) only defines a conserved 
quantity for a Killing vector , the vector fields (a e {0, 1,2,3}) corre- 
sponding to the 4-translations at get closer to an infinitesimal symmetry 
as one approaches Therefore, one can expect that a Komar- like expres- 
sion makes sense for a given cross-section Su of This is indeed the case 
and the evaluation of the integral does not depend on how we get to Su- 



(48) 
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However, the integral does depend on the representative in the class of 

vectors corresponding to the asymptotic symmetry. This is cured by imposing 
a divergence- free condition on [40] . Bondi-Sachs 4-momentum at Su C 
is then defined as 



Alternatively, ambiguities in the Komar integral can be solved by dropping 
the condition on the divergence and adding a term ccV^^^ to the surface 
integral. When a = 1 the resulting integral is called the linkage [41]. The 
discussion of Bondi-Sachs angular momentum is more delicate. We refer the 
reader to the discussion in section 3.2.4 of [1]. 

Bondi energy and positivity of gravitational radiation energy 

Bondi energy Ess (the zero component of the Bondi-Sachs 4-momcntum) is 
a decreasing function of the retarded time. More concretely, Bondi energy 
satisfies a loss equation 



where F > can be expressed in terms of the squares of the so-called news 
functions. In [21] it is shown that, if the news tensor satisfies the appropriate 
conditions, then Bondi mass coincides initially with the ADM mass (see also 
[22]). Bondi energy is interpreted as the remaining of the ADM energy in 
the process of energy extraction by gravitational radiation. As for the ADM 
mass, a positivity result holds for the Bondi mass [42, 43]. These properties 
constitute the underlying conceptual/structural justification of our under- 
standing of energy radiation by gravitational waves: gravitational radiation 
carries positive energy away from isolated radiating systems, and the total 
radiated energy cannot be bigger than the original total ADM energy. 

4 Notions of mass for bounded regions: quasi-local 
masses 

As commented in section 1.2, the convenience of associating energy-momentum 

with the gravitational field in given regions of the spacctimc is manifest in 
very different contexts of gravity physics. More specifically, mathematical 
and numerical General Relativity or approaches to Quantum Gravity pro- 
vide examples where we need to associate such an energy-momentum with a 
finite region of spacetime. This can be either motivated by the need to de- 
fine appropriate physical/astrophysical quantities, or by the convenience of 
finding quasi-local quantities with certain desirable mathematical properties 
(e.g. positivity, monotonicity...) in the study of a specific problem. 





(49) 




(50) 
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There exist many different approaches for introducing quasi-local pre- 
scriptions for the mass and angular momentum. Some of them can be seen 
as quasi-localizations of successful notions for the physical parameters of the 
total system, such as the ADM mass, whereas other attempts constitute gen- 
uine ab initio methodological constructions, mainly based on Lagrangian or 
Hamiltonian approaches. An important drawback of most of them in the con- 
text of the present article is that, typically, they involve constructions that 
are difficult to capture in short mathematical definitions without losing the 
underlying physical/geometrical insights. An excellent and comprehensive re- 
view is reference [1] by Szabados. 

Ingredients in the quasi-local constructions. 

First, the relevant bounded spacetime domain must be identified. Typically, 
these are compact space-like domains D with a boundary given by a closed 
2-surface S. Explicit expressions, such as relevant associated integrals, are 
formulated in terms of either the (3-dimensional) domain D itself or on its 
boundary S. In particular, conserved-current strategies permit to pass from 
the 3-volume integral to a conserved-chaxge-like 2-surface integral. In other 
cases, 2-surface integrals are a consequence of the need of including boundary 
terms for having a correct variational formulation (as it was the case in the 
Hamiltonian formulation of section 3.3). 

We have already presented an example of quasi-local quantity in sec- 
tion 2, namely the Komar quantities. Since symmetries will be absent in the 
generic case, an important ingredient in most quasi-local constructions is the 
prescription of some vector field that plays the role that infinitesimal symme- 
tries had played in case of being present. In connection with this, one usually 
needs to introduce some background structure that can be interpreted as a 
kind of gauge choice. 

Finally, different plausibility criteria for the assessment of the proposed 
quasi-local expressions (e.g. positivity, monotonicity, recovery of known lim- 
its...) need to be considered (see [1]). 

4.1 Some relevant quasi- local masses 
Round spheres. Misner-Sharp energy 

In some special situations, as it is the case of isolated systems above and 
some exact solutions, there is agreement on the form of the gravitational field 
energy-momentum. Another interesting case is that of spherically symmetric 
spacetimes, where the rotation group SO{3) acts transitively as an isometry. 
Orbits under this rotation group are round spheres <S. Then, using the areal 
radius as a coordinate (47rr^ = ^): an appropriate notion of mass/energy 
was given by Misner and Sharp [44] 
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E{S) := -r\R^, 



2 /J,!/ 2 per 



(51) 



where ^e^,y = nf s^^pafiv (cf. section 1.3) is the volume element on S. This 
expression is related to the so-called Kodama vector K'^, that can be defined 
in spherically symmetric spacetimes and such that W ^{C^" K^,) = 0. The 
current = G^'^K^ is thus conserved and, taking D as a solid ball of radius 
Ta , the flux of through the round boundary dD actually equals the change 
in time of the mass expresion (51). Misner-Sharp proposal is considered as 
the standard form of quasi-local mass for round spheres. 

Brown- York energy 

The rationale of the approach in Ref. [45] to quasi-local energy strongly 
relies on the well-posedness of a variational problem for the gravitational 
action. The adopted variational formulation is essentially the one outlined in 
section 3.3 (where the discussion was in fact based in the treatment in [3] 
adapted from [45]). However, if the main interest is placed in the expressions 
of quasi-local parameters and not in the details of the symplectic geometry 
of the system phase space, a full Hamiltonian analysis does not need to be 
undertaken and one can rather follow a Hamilton-Jacobi one. The latter 
starts from action (21) defined on the spacetime domain V. We recall that 
the boundary dV is given by two space-like hypersurfaces Si and S2 and 
a time-like tube B, such that the 2-spheres Si are the intersections between 
Si and B. The metric and extrinsic curvatures on Si are given by 7^^ and 
Rf^"^, whereas those on B are denoted by and Pf^" . A Hamilton-Jacobi 
principal function can then be introduced by evaluating the action S on 
classical trajectories. An arbitrary function S*" of the data on the boundaries 
can be added to S [it is the responsible of the reference terms with subindex 
in expression (21)]. The principal function is given by Sc\ ■= {S — S^) jci 
and Hamilton-Jacobi equations are obtained from its variation with respect 
to the data at the final slice S2. One of the Hamilton-Jacobi equations leads 
to the definition of a surface stress- energy-momentum tensor as 



This tensor satisfies a conservation-like equation with a source given in terms 
of the matter energy- momentum tensor T'*'^. This motivates the definition of 
the charge Qs{^^) associated with a vector as 




2_^ 




(52) 




(53) 



whose change along the tube B is given by a matter flux. This expression is 
analogous to (1) in the matter case (here S C B and is the time-like unit 
normal to S and tangent to B). 
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Using the 2+1 decomposition induced by a 3+1 space-like slicing {St}j 
wc can decompose the tensor r^" as we did for the matter energy-momentum 
tensor T^" in Eq. (12). Writing exphcitly the time-hke components, it results 



1 

Stt 



Expressing the vector in the 3+1 decomposition = ^n^ + and con- 
sidering a 2-surface S lying in a slice of {St}, we have 

Qsin = i ipT'^n^^q Sx =1 {ie- ^Jp) ^ cfx . (55) 
Js Js 

The Brown- York energy is then [cf. with the ADM mass expression (41)] 

EnY{S,n^):=Qs{n^) = -^jj,H-Ho)^dPx . (56) 

Note that this expression exphcitly depends on the manner in which <S is 

inserted in some space-like 3-sliee. In this sense, it corresponds to an energy 
(depending on a boost) rather than a mass. 

Kijowski, Epp, Liu-Yau and Kijowski-Liu-Yau expressions 

We briefly comment on some expressions that can be related to the Brown- 
York energy. Studying more general boundary conditions than the ones in 
[45], Kijowski proposed the following quasi-local expression for the mass [46] 

^KiJ := / ^^^^^^^-^^ d^x , (57) 



167r Js Ho 

where H = iJ^^g^"^ and L = L/j^q'^" are the traces of the extrinsic curvatures 
of <S with respect to unit orthogonal space-like s'^ and time-like n'^ vectors, 
i.e. n^Sfj, = (cf. notation in section 1.3). Apart from the choice of the 
background terms Hq, this expression only depends on S, and not in the 
manner of embedding it into some space-like hypersurface. Using a different 
set of boundary conditions, another quasi-local quantity was introduced by 
Kijowski (referred to as a free energy). The same quantity was later derived 
by Liu and Yau, using a different approach [47] . We will refer to the resulting 
quasi-local energy as the Kijowski-Liu-Yau energy, having the form 



(58) 



On the other hand, aiming at removing the dependence of Brown- York en- 
ergy on the space-like hypersurface, Epp [48] proposed the following boost- 
invariant expression 
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:= ^ £ {Vm^ - {L^f - - L^) ■ (59) 

Note that Brown- York energy can be seen as a gravitational field version 
of the quasi- local matter energy (6), whereas Epp's expression rather corre- 
sponds to the matter mass (7). For further recent work along this approach 
to quasi-local mass, see [49, 50]. 

Hawking, Geroch and Hayward energies 

Hawking energy. 

Given a topological sphere 5, its Hawking energy is defined as [51] 

(£) (k) 

where 9^ = q^'^Ofiu and 9- = q^'^OfiJ are the expansions associated with 
outgoing and ingoing null normals (cf. notation in section 1.3). It can be 
motivated by understanding the mass surrounded by the 2-sphere S as an 
estimate of the bending of ingoing at outgoing light rays from S. An average, 
boost-independent measure of this convergence-divergence behaviour of light 
rays is given by §g6+6-\. Then, from the Ansatz A + B §g9+9-\, the 
constants A and B are fixed from round spheres in Minkowski and from the 
horizon sections in Schwarzschild spacetime. 

Hawking energy depends only on the surface S and not on any particular 
embedding of it in a space-like hypersurface. In the spherically symmetric case 
it recovers the standard Misner-Sharp energy (51). For apparent horizons, or 
more generally for marginally trapped surfaces, it reduces to the irreducible 
mass accounting for the energy that cannot be extracted from a black hole 
by a Penrose process and that is given entirely in terms of the area. Hawking 
energy does not satisfy a positivity criterion, since it can be negative even 
in Minkowski spacetime. However, for large spheres approaching null infinity, 
E^{S) recovers Bondi-Sachs energy, whereas for spheres approaching spatial 
infinity it tends to the ADM energy. Though it is not monotonic in the 
generic case, monotonicity can be proved for sequences of spheres obtained 
from appropriate geometric flows. This has a direct interest for the extension 
of Huisken & Ilmanen proof [52] of the Riemaniann Penrose inequality to the 
general case. 

Geroch energy. 

For a surface <S embedded in a space-like hypersurface S, Geroch energy [53] 
is defined as 
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where H is again the trace of the extrinsic curvature of S inside E. Geroch 
energy is never larger than Hawking energy, but it can be proved that it also 
tends to the ADM mass for spheres approaching spatial infinity. 

The relevance of Geroch energy lies on its key role in the first proof of 
the Riemaniann Penrose, by Huiskcn & Ilmancn [52] (see also section 5.1). In 
particular, use is made of the monotonicity properties of Eq under an inverse 
mean curvature flow in S. 

Hayward energy. 

Some generalizations of Hawking energy exist. A vanishing expression for flat 
spacetimes can be obtained by considering the modified expression 

where the shears cr^^ and cr"^ are the traceless parts of 6*^^^ and O^^J , respec- 
tively, still asymptotes to the ADM energy at spatial infinity, but does not 
recover Bondi-Sachs energy at null infinity (but rather Newman-Unti one; cf. 
references in [1]). Related to this modified Hawking energy, Hayward has pro- 
posed [54] another quasi-local energy expression by taking into account the 
anholonomicity form i?^^ , one of the normal fundamental 1-forms introduced 
in section 1.3 

Though the divergence-free part of i?^ can be related to angular momentum 
(see below), this 1-form is a gauge dependent object changing by a total 
differential under a boost transformation. Therefore, some natural gauge for 
fixing the boost freedom is needed. 

Bartnik mass 

Bartnik quasi-local mass is an example of quasi-localization of a global quan- 
tity, in particular the ADM mass. In very rough terms, the idea in Bartnik's 
construction consists in defining the mass of a compact space-like 3-domain 
D as the ADM mass of that asymptotically Euclidean slice S that contains 
D without any other source of energy. The strategy to address this absence 
of further energy is to consider all plausible extensions of D into Euclidean 
slices, calculate the ADM mass for all them, and then consider the infimum 
of this set of ADM masses. In more precise terms, let us consider a compact, 
connected 3-hypersurface D in spacetime, with boundary S and induced met- 
ric jij. Bartnik's construction actually focuses on time-symmetric Kij = 
domains D. Let us also assume that a dominant energy condition (though 
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the original formulation in [55] makes use of a weak-energy-constraint con- 
dition) is satisfied. In a time-symmetric context this amounts to the posi- 
tivity of the Ricci scalar, "^i? > 0. One can then define P{D) as the set of 
Euclidean time-symmetric initial data sets {U.^-jij) satisfying the dominant 
energy condition, with a single asymptotic end, finite ADM mass Madm{S), 
not containing horizons (minimal surfaces in this context) and extending D 
through its boundary S. Then, Bartnik's mass [55] is defined as 



The no-horizon condition is needed to avoid extensions with arbi- 

trarily small ADM mass. There is also a spacetime version of Bartnik's con- 
struction, not relying on an initial data set on D but only on the geometry of 
2-surfaces S. Let us define 'P{S) as the set of globally hyperbolic spacetimes 
{M.,g^v) satisfying the dominant energy condition, admitting an asymptoti- 
cally Euclidean Cauchy hypcrsurfacc E with finite ADM mass, not presenting 
an event horizon and such that S is embedded (i.e. both its intrinsic and ex- 
trinsic geometry) in {M.,g^v)- Then, one defines 



The comparison between Mb(-D) and M^idD) is not straightforward, due 
to issues regarding the horizon characterisation. Prom the positivity of the 
ADM mass it follows the non- negativity of the Bartnik mass Mb{D). In fact, 
Mb(-D) = characterises D as locally flat. From the definition (64) it also 
follows the monotonicity of MB(r>), i.e. if Di C D2 then MB(r>i) < Mb{D2). 
Bartnik mass tends to the ADM mass, as domains D tend to Euclidean slices 
(the proof makes use of the Hawking energy introduced above). Another 
interesting feature, consequence of the proof of the Riemannian Penrose con- 
jecture [52], is that Bartnik mass reduces to the standard form E{S) in Eq. 
(51) for round spheres. However, the explicit calculation of the Bartnik mass 
is problematic. An approach to its practical computability is provided by 
Bartnik's conjecture stating that the infimum in (64) is actually a minimum 
realised by an element in 'P{D) characterised by its stationarity outside D. 
Further developments of these ideas have been proposed by Bray (cf. [56]). 

Before concluding this subsection, we mention the explicit construction 
in [57] of quantum analogs for some of the previous quasi-local gravitational 
energies (specifically for Brown- York, Liu-Yau, Hawking and Geroch ener- 
gies), as operators acting on the appropriate representation Hilbert space in 
a particular approach to quantum gravity (namely, loop quantum gravity). 

4.2 Some remarks on quasi- local angular momentum 

Spinorial techniques provide a natural setting for the discussion of angular 

momentum. This does not only apply to angular momentum, since spino- 
rial and also twistor techniques define a framework where further quasi-local 



MYi{D) inf {Madm(^), such that E e V{D)} 



(64) 



Mb(5) := inf {Madm(A^), such that M e V{S)} 



(65) 
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mass notions can be introduced (e.g. Penrose mass), and known results can 
be reformulated in particularly powerful formulations (e.g. the discussion of 
positive mass theorems using the Nester-Witten form). However, in this ar- 
ticle we will not discuss these approaches and we refer the reader to the 
relevant sections in Ref. [1]. We will focus on certain aspects of quasi-local 
expressions for angular momentum of Komar-like type. As it was shown in 
section 2.2, choosing a two-sphere <S in a 3-slice adapted to the axial symme- 
try (/)^, a 1-form can be found such that the Komar angular momentum 
is expressed as 



In particular, in the Komar expression (18) we have = Qf^^KiypsP, whereas 
in the spatial infinity expression (45) this is modified by a term proportional 
to the trace K of the extrinsic curvature. The same applies for an angu- 
lar momentum defined from the Brown- York charge (53) when plugging the 
expression for in (54) into (55), where (j)'^ does not need to be a sym- 
metry. The normal fundamental 1-forms Ofi on «S (cf. section 1.3) provide 
another avenue to L^. In this section we assume the form (66) for the an- 
gular momentum and comment on some approaches to the determination of 
the {quasi- symmetry) axial vector cj)!^. 

Divergence-free and quasi-Killing axial vectors 

No ambiguity for (j)^^ is present when an axial symmetry exists on S: (j)^ is 
taken as the corresponding Killing vector. In the absence of such a symmetry, 
we must address two issues. First, expression (66) depends on the space-like 
3-slice in which S is embedded. This follows from the modification of the 1- 
form by a total differential under a boost transformation: L^-^-'^Dnf 
(cf. boost/normalization transformation of f?^^^ in section 1.3). Angular mo- 
mentum can be associated with <S, independently of any hypersurface S, 
by demanding the axial vector to be divergence- free: '^D^cj)'^ = 0. Then, the 
boost-induced modification vanishes under integration. Second, the physical 
meaning of J{(t)^ ) is unclear if references to a symmetry notion are completely 
dropped. In this sense, different approaches exist aiming at defining appropri- 
ate quasi-Killing notions. We simply mention here some recent works along 
these lines. In the context of isolated horizons (see next subsection) a pre- 
scription for the determination of a quasi-Killing axial vector on black hole 
horizons has been proposed in [58], though the divergence-free character is 
not guaranteed. Ref. [59] presents an approach for finding an approximate 
Killing vector by means of a minimization variational prescription that re- 
spects the divergence- free character of </>'*. In the context of dynamical or 
trapping horizons [60, 61, 62], a unique divergence-free vector </>'* can be cho- 
sen such that it is preserved by the unique slicing of the (space-like) horizon 
worldtube by marginally outer trapped surfaces [63]. Also in the context of 




(66) 
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dynamical horizons, a proposal for (j)^ has been made in [64] relying on a con- 
formal decomposition of the metric on S. Sec Rcf. [65] for a discussion 
of the divergence-free character of vector fields associated with quasi-local 
observables on S. 

Equation (66) only provides the expression for the component of the an- 
gular momentum vector that is associated with the vector 0^. If we are inter- 
ested in determining the total angular momentum vector, a sensible prescrip- 
tion for the other two components is needed. This is an important practical 
issue in numerical simulations (see e.g. [66]). 

4.3 A study case: quasi-local mass of black hole isolated horizons 

The need of introducing some additional structure has been discussed above 
in diflFerent settings (e.g. symmetries for Komar quantities and asymptotic 
flatness for ADM and Bondi masses) . Wc illustrate now this issue in a quasi- 
local context related to equilibrium black hole horizons. 

A brief review of isolated horizons 

The isolated horizon framework introduced by Ashtekar and collaborators 

[62] provides a quasi-local setting for characterising black hole horizons in 
quasi-equilibrium inside an otherwise dynamical spacetime. It presents a hi- 
erarchical structure with different quasi-equilibrium levels. The minimal no- 
tion of quasi-equilibrium is provided by the so-called non- expanding horizons 
(NEH). Given a Lorentzian manifold, a NEH is a hypcrsurface T-L such that: 

i) "H is a null hypersurface of topology 5^ x M that is sliced by marginally 
(outer) trapped surfaces, i.e. the expansion of the null congruence associ- 
ated with the null generator vanishes on %: 6^^^ = q'^^O^^l = 0. 

ii) Einstein equation is satisfied on H. 

iii) The vector —T^^^f is future directed. 

The geometry of a NEH is characterised by the pair (g^y, V^), where q^v is 
the induced null metric on % and is the unique connection (not a Levi- 
Civita one) induced from the ambient spacetime connection. characterises 
the extrinsic geometry of the NEH. A certain combination of components in 
can be put together to define an intrinsic object on %, namely the 1-form 
Wfj, characterised by 

V^r=w^t . (67) 

Defining a surface gravity as := i'^to^, the acceleration expression for 

is given by: Vgl^ = K,(^t)^^- On the other hand, the projection of 
on S recovers the fundamental normal 1-form: J?^^^ = q^''u)p. The quasi- 
equilibrium hierarchy is introduced by demanding the invariance of the null 
hypersurface geometry under the (evolution) fiow in a progressive manner: 
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1. A NEH is characterised by the time-invariance of the intrinsic geometry 

2. A weakly isolated horizon (WIH) is a NEH, together with an equiva- 
lence class of null normals [^^], for which the 1-form ui^ is time-invariant: 
jC^cOh = 0. This is equivalent to the (time and angular) constancy of the 
surface gravity: V^K(^) — 0. 

3. An isolated horizon (IH) is a WIH on which the whole extrinsic geometry 
is time-invariant: [£|^, V^] = 0. 

The NEH and IH quasi-cquilibrium levels represent genuine restrictions on 
the geometry of 'H as a hypersurface in the ambient spacetime. On the con- 
trary, a WIH structure can always be implemented on a NEH by an appro- 
priate choice of the null normal £^ normalization. In this sense, a WIH does 
not represent a higher level of quasi-equilibrium than a NEH. However, from 
the point of view of the Hamiltonian analysis of spacetimes with a black 
hole in quasi-equilibrium as an inner boundary, the WIH notion proves to 
be crucial for the correct definition of the phase space symplectic structure 
and, more concretely, for the sound formulation of the quasi-local mass and 
angular momentum of the horizon. 

An overview of the Hamiltonian analysis of isolated horizons 

Conserved quantities under horizon symmetries 

As in the presentation of ADM quantities in section 3.3, mass and angular 
momentum of isolated horizons are introduced as conserved quantities under 
appropriate symmetries (see [67, 68] and the outline in Appendix C of [69] for 
further details on the following discussion). One starts from a symmetry of the 
horizon structure in the Lorentzian spacetime manifold and then constructs 
an associated canonical transformation in the phase or solution space of the 
system. The conserved quantity under this canonical transformation provides 
the relevant physical quantity. In view of the variational problem (see below), 
a WIH is the relevant horizon structure to be considered in this context. A 
vector field preserves the WIH structure (W is a WIH-symmetry) if 

r^^** = const • , C^q^^ = , = . (68) 

WIH-symmetries are of the form = cw^'^ + bwS'^, where cw and bw are 
constants on H and is a Killing vector of any spatial section <S of H. 

Variational problem for spacetimes containing WIHs 

In order to set up the Hamiltonian treatment, we need first to define a well- 
posed variational problem. Here we are interested in the variational problem 
for asymptotically flat spacetimes containing a WIH. We will furthermore 
demand this WIH to contain an axial Killing vector (p^. The variational 
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problem is then set in the region contained between two asymptotically Eu- 
clidean slices i7_ and E^, spatial infinity i'^ and the part of TL between an 
initial horizon slice S- = T-LD and a final one Sj^ =T-LC\ The action, 
as in Eq. (21), can be written [67, 68] as the sum of a bulk and a bound- 
ary term at spatial infinity, and the variation of the dynamical fields is set 
to vanish on the slices S- and S+. No boundary term associated with the 
inner boundary % is introduced. The variational problem is well-posed, in 
particular the Einstein equation is recovered, as long as the condition 



holds, where uj^ has been introduced in (67) and is the vohmic 2-form on 
sections S ofH. The crucial ingredient in the well-posedness of the problem 
is precisely the WIH structure. This is the additional structure needed in 
order to guarantee the vanishing of (69), so that the variational problem is 
correctly posed and quasi-local quantities can be defined. 

Phase space, canonical transformations and physical quantities 

The phase space is defined by the couple (-T, J) where F is an infinite- 
dimensional manifold where each point represents a solution to the Einstein 
equation containing a WIH, and J is a symplectic form (a closed 2-form) on 
r in terms of which the Poisson bracket is defined. In particular, a vector 
field X on r generates a canonical transformation if it leaves the symplectic 
form invariant: Cx"^ — 0- Using the closedness of J this is equivalent to the 
exactness of the 1-form ixJ, i-e. to the (local) existence of a function Hx 
such that ixJ = SHx (where S denotes the differential in F). In particular, 
the quantity Hx defined on the phase space is preserved along the flow of X. 
In this context, first, the symplectic form can be obtained from the ac;tion by 
using the conserved symplectic current method [70] and, second, a vector field 
X\Y on F can be constructed from a WIH-symmetry on H (cf. [67, 68] 
for details) . For the correct definition of a physical parameter associated with 
a given WIH-symmetry W^, we must assess if the corresponding Xw pre- 
serves the canonical form, i.e. if J is locally exact. If this is the case, the 
conserved quantity is simply read from the associated explicit expression of 
the Hamiltonian • When this scheme is applied to the axial symmetry 
(j)'^ onH, the corresponding X^ turns out to be automatically an infinitesimal 
canonical transformation and the conserved quantity has the form 



where St is any spatial section of H. This prescription for J-}i exactly coincides 
with the Komar expression (18). The mass discussion is more subtle. In this 
case the WIH-symmetry t^ associated with time evolution is chosen as an 




(69) 




(70) 
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appropriate linear combination of the null normal £^ and the axial vector (^'*. 
It is then found 

ix,3 = SEadm - SAh + ^2(4) SJn) , (71) 

where K(^t) and i7(t) are functions on F and A-u and J-u correspond, respec- 
tively, to the area of any section of H and to the horizon angular momentum 
in (70). The right-hand-side expression is (locally) exact if functions K(t) and 

depend only on A-u and J-^, and satisfy: = Sir ^J^^ ■ A function 
only depending on A-^ and J-h then exists, such that we can write 

6El^ = ^(*)(^^' -^^^ SAn + 0^t){An, Jn) SJn ■ (72) 

To finally determine the quasi- local mass M^, the functional form of E\^{A-u, J-u) 
is normalized to the one in the stationary Kerr family in F. Note that this is 
only justified once E^ has been shown to depend only on A-u and J-u , a non- 
trivial result. In sum, for isolated horizons M-u,{A-u,, Jn) '■— AfKerr(^w, •/«), 
given by the Christodoulou mass expression [71] . 

Remark 8 (Quasi-local first law of black hole dynamics). Expression (72) ex- 
tends the first law of black hole dynamics (see section 5.2) from the stationary 
setting to dynamical spacetimes where only the black hole horizon is in equi- 
librium. 



5 Global and quasi-local quantities in black hole physics 

As an application, we briefly comment on some relevant issues concerning 
mass and angular momentum in the particular case of black hole spacetimes. 

5.1 Penrose inequality: a claim for an improved mass positivity 
result for black holes 

In the context of the established gravitational collapse picture, Penrose [16] 
proposed an inequality providing an upper bound for the area of the spatial 
sections of black hole event horizons in terms of the square of the ADM mass. 
This conjecture followed from a heuristic chain of arguments including rigor- 
ous results (singularity and black hole uniqueness theorems), together with 
conjectures such as weak cosmic censorship and the stationarity of the final 
state of the evolution of a black hole spacetime. A local-in-time version of the 
Penrose inequality can be formulated in terms of data on a Euclidean slice. 
In this version Penrose conjecture states that, given an asymptotically Eu- 
clidean slice S containing a black hole under the dominant energy condition, 
the following inequality should be satisfied 
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^min < leTTMloM , (73) 

where A^i^ is the minimal area enclosing the apparent horizon. In addition, 
equality is only attained by a slice of Schwarzschild spacetime. Though this 
was originally proposed in an attempt to construct counter-examples to the 
weak cosmic censorship conjecture, growing evidence has accumulated sup- 
porting its generic validity. Beyond spherical symmetry [72], a formal proof 
only exists in the Riemannian case, JC^^ = 0, where the original derivation 
[52] (see also [73]) makes use of some of the quasi-local expressions presented 
in section 4 (cf. discussion about Geroch and Hawking energies, that coincide 
in this time-symmetric case K^^i, = 0). The intrinsic geometric relevance of 
the Penrose inequality is reflected in its alternative name as the isoperimetric 
inequality for black holes [74] . 

Penrose inequality can also be seen as strengthening the positive ADM 
mass theorem in section 3.3, for the case of black hole spacetimes: the ADM 
mass is not only positive but must be larger than a certain positive-definite 
quantity. Though it is tempting to identify this positive quantity with some 
quasi-local mass associated with the black hole, e.g. with its irreducible mass 
A =: 16ttM?j. related to the Hawking mass (60), a caveat follows from the 
fact that the relevant minimal surface of area j4,nin does not necessarily co- 
incide with the apparent horizon, as examples in [75] show. In any case, this 
geometric inequality represents a bridge between global and quasi-local prop- 
erties in black hole spacetimes and has become one of the current geometric 
and physical/conceptual main challenges in General Relativity. 

5.2 Black hole (thermo-) dynamics 

A set of four laws was established in [76] for stationary black holes. These 

black hole laws are analogous in form to the standard thermodynamical laws. 
Though this analogy is compelling, the fundamental nature of such relation 
was only acknowledged under the light of Hawking's discovery [77] of the 
(scmiclassical) thermal emission of particles from the event horizon (Hawking 
radiation). Given a stationary black hole spacetime with stationary Killing 
vector t^, black hole rigidity theorems [78] imply the existence of a second 
Killing vector that coincides with the null generators on the horizon. 
We can write k^ =t^ + Qh<P^ , where (f)^ is an axial Killing vector and Qh is 
a constant referred to as the angular velocity of the horizon (see also [4]). We 
can write k'^'Vi^k^ = Kk^ on the horizon, which defines the surface gravity 
function k. The zeroth law of black hole mechanics then states the constancy 
of the surface gravity on the event horizon. The second law, namely Hawking's 
area theorem [79, 80], guarantees that the area of the event horizon never 
decreases, whereas the third law states that the surface gravity k, cannot be 
reduced to zero in a finite (advanced) time (see [81] for a precise statement). 
In the present context, we are particularly interested in the first law, since 
it relates the variations of some of the quasi-local and global quantities we 
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have introduced in the text, in the particular black hole context. First law 
provides an expression for the change of the total mass M of the black hole (a 
well-defined notion since we deal with asymptotically flat spacetimes) under 
a small stationary and axisymmetric change in the solution space 



where A is the area of a spatial section of the horizon, and Jh is the Komar 
angular momentum associated with the axial Killing t/)^. Equation (74) re- 
lates the variation of a global quantity M = Madm at spatial infinity on the 
left- hand- side, to the variation of quantities locally defined at the horizon, 
on the right-hand-side. In particular, we could express the variation of the 
horizon area in terms of the variation of the irreducible local mass Mi„, as 
SA = 327rMirr(5Mirr. Such a formulation actually plays a role in the critcirion 
for constructing sequences of binary black hole initial data corresponding to 
quasi-circular adiabatic inspirals (cf. [30] and the first law of binary black 
holes in [82]). Derivation of (74) involves the notions of ADM mass, as well 
as the generalization to stationarity of the Smarr formula for Kerr mass [stat- 
ing M = 2QhJh + KA/(47r)] by using the Komar mass expression. Result 
(72) in section 4.3 provides an extension of this law to black hole spacetimes 
non-necessarily stationary, but containing an isolated horizon for which an 
unambiguous notion of black hole mass can be introduced. Quasi-local at- 
tempts to extend the first law to the fully dynamical case have been explored 
in the dynamical and trapping horizon framework [60, 61, 62, 83]. However, 
the lack of a general unambiguous notion of quasi-local mass prevents the 
derivation of a result analogous to (74) or (72), i.e. the equality between 
the variation of two independent well-defined quantities. In the quasi-local 
dynamical context, an unambiguous law for the area evolution can be deter- 
mined (see e.g. [84, 85] and references therein). The latter can then be used 
to define a flux of energy through the horizon by comparison with (74). 

5.3 Black hole extremality: a mass-angular momentum inequality 

Subextremal Kerr black holes are characterised by presenting angular mo- 
menta bounded by their total masses. Keeping axisymmetry, it has been 
recently shown [86, 87, 88] (see also [89]) that the inequality 



holds also for vacuum, maximal (K=0), axisymmetric Euclidean data. More- 
over, equality only holds for slices of extremal Kerr. Inequality (75) provides 
a non-trivial relation for black hole spacetimes between precisely the two 
physical quantities we are focusing on in this review. It is natural to explore 
if some analogous inequality holds when moving away from axisymmetry and 
when considering only the local region around the black hole. Attempts have 




(74) 



JkI < Ml 



ADM 5 



(75) 
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been done in this sense, but they all must face the ambiguities resulting from 
the absence of canonical expressions for quasi-local masses and angular mo- 
menta. In order to illustrate the caveats to keep in mind when undertaking 
this kind of discussion, one can consider the case in which Komar quantities 
arc used for constructing a truly quasi-local analogue of expression (75) for 
axisymmetric stationary data: initial data have been constructed [90] where 
the quotient | Jk|/M^ on the black hole horizon can become arbitrarily large. 
Interestingly, these studies have led to the formulation [91] in axisymmetry of 
the related conjecture 87r| Jk| < A, only involving intrinsic quantities on the 
horizon. This inequality has been proved to hold in the stationary axisym- 
metric case [92], as well as in a generalization including the electromagnetic 
field and the associated electric charge on the horizon [93]. 

6 Conclusions 

The problem of characterising the energy-momentum and angular momentum 
associated with the gravitational field in General Relativity has been present 
since the birth of the theory and controversies have plagued its already long- 
standing history. Understanding that no local density of energy-momentum 
can be identified for the gravitational field has challenged the validity of 
the mass and angular momentum cherished notions from non-gravitational 
physics, when trying to perform a straightforward extension of these concepts 
to the gravitational field in a general relativistic setting. 

The study of specific problems suggests concrete and/or partial solutions. 
In this spirit, at low velocities and weak self-gravities post-Newtonian ap- 
proaches handle consistent notions of mass and angular momentum and the 
same holds in perturbative approaches around exact solutions, for which 
physical quantities can be identified unambiguously. In the same line, a 
(quasi-local) notion of the energy carried by a gravitational wave can be 
introduced as an average along the wavelength, proving to provide a use- 
ful notion in practical applications. A particular setting of singular concep- 
tual importance is that of isolated systems in General Relativity, specifically 
through their characterisation as asymptotically flat spacetimes. The notions 
of total ADM and Bondi-Sachs energy-momentum provide well-defined quan- 
tities that, on the one hand, have clarified important conceptual issues such 
as the capability of gravitational waves to actually carry energy away from 
a system and, on the other hand, they also represent inestimable tools in 
practical applications due to their intrinsic/geometric character. Positivity 
theorems for the total mass represent without any doubt some of the most 
important and profound resiilts in General Relativity. The combination of the 
success in isolated systems, together with the absence of a gravitational local 
energy-momentum density, has led to the consideration that the whole effort 
for the search of a local expression for the gravitational energy represents an 
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ill-posed or pseudo-problem (see e.g. [11]). But at the same time, and mo- 
tivated by practical needs and/or fundamental physical reasons (cf. in this 
sense [94] for a related discussion on quasi-local issues regarding observables 
in Quantum Field Theory), important efforts have been devoted to the intro- 
duction of quasi-local notions of gravitational energy-momentum associated 
with extended but finite regions of the spacetime. In this respect, significant 
insights into the structure of the gravitational field have been achieved, with 
applications in diverse conceptual and practical contexts. But it must be ac- 
knowledged, as it is referred in [1], that the status of the quasi-local mass 
studies is in a kind of post-modem situation in which the devoted intensive 
efforts have resulted in a plethora of proposals with no obvious definitive and 
entirely satisfying candidate. 

A moderate (intermediate) position that avoids radical skepticism against 
the quasi-local approach would consist in assuming that mass and (in a more 
restricted sense) angular momentum can be unambiguously defined only as 
global quantities for isolated systems. But accepting, at the same time, that 
quasi-local expressions provide meaningful and insightful quantities that are 
inextricably subject to the need of making systematically explicit the specific 
setting in which they are defined (one can make the analogy with the notion 
of effective mass in solid state physics, where different masses can be simulta- 
neously employed for the same particle as long as their specific purposes are 
clearly stated^). The moral of the whole discussion in this article is that the 
formulation of meaningful global or quasi-local mass and angular momentum 
notions in General Relativity always needs the introduction of some addi- 
tional structure in the form of symmetries, quasi-symmetries or some other 
background structure. This point must be explicitly kept in mind whenever 
employing the so-defined physical quantities, specially when extrapolating or 
performing compared analysis. 

Acknowledgements. The authors wish to thank the organisers of the Or- 
leans School on Mass for their kind invitation and encouragement. 



References 

1. Szabados, L.B.: Quasi-local energy-momentum and angular momentum in gr: 
a review article. Liv. Rev. Rclat. 7, 4 (2004). URL (cited on 8 May 2009): 
http://www.livingreviews.org/lrr-2004-4 

2. Gourgoulhon, E.: 3-1-1 formalism and bases of numerical relativity, lectures 
delivered at Institut Henri Poincare in 2006, available as gr-qc/0703035 

3. Poisson, E.: A Relativist's Toolkit, The Mathematics of Black-Hole Mechanics. 
Cambridge University Press, Cambridge (2004) 

4. Wald, R.M.: General Relativity. Chicago University Press (1984) 



^We thank B. Carter for his many comments and insights in this discussion, 
and in particular for bringing us to the solid state analogy for gravitational masses. 



Mass and Angular Momentum in General Relativity 



37 



5. Senovilla, J.M.M.: Super-energy tensors. Class. Quant. Grav. 17, 2799-2842 

(2000). do: 10.1088/0264-9381/17/14/313 

6. Noether, E.: Nachr. Gcs. Wiss. Goettingen 2, 235 (1918) 

7. Bclinfantc, J.: Physica 6, 887 (1939) 

8. Belinfante, J.: Physica 7, 449 (1940) 

9. Greiner, W., Reinhardt, J.: Field Quantization. Springer- Verlag, Berlin (1996) 

10. Chang, C.C., Nester, J.M., Chen, CM.: Pseudotensors and quasilocal 
energy-momentum. Phys. Rev. Lett. 83(10), 1897-1901 (1999). DOI 
10.1103/PhysRevLett.83.1897 

11. Misner, C.W., Thorne, K.S., Wheeler, J.A.: Gravitation. Freeman (1973) 

12. Wald, R.M.: The thermodynamics of black holes. Living Reviews in Relativity 
4(6) (2001). URL http://www.livingrcvicws.org/hr-2001-6 

13. Jaramillo, J.L., Valiente Kroon, J. A., Gourgoulhon, E.: From Geometry 
to Numerics: interdisciplinary aspects in mathematical and numerical rel- 
ativity. Class. Quantum Grav. 25, 093,001 (2008). DOI 10.1088/0264- 
9381/25/9/093001 

14. Komar, A.: Covariant conservation laws in general relativity. Phys. Rev. 113, 

934 (1959) 

15. Katz, J.: A note on Komar's anomalous factor. Class. Quantum Grav. 2, 423 

(1985) 

16. Penrose, R.: Naked singularities. Annals N. Y. Acad. Sci. 224, 125 (1973) 

17. Penrose, R.: Zero rest-mass fields including gravitation: asymptotic behaviour. 
Proc. Roy. Soc. Lond. A 284, 159 (1965) 

18. Geroch, R.P.: Structure of the gravitational field at spatial infinity. J. Math. 
Phys 13, 956-968 (1972) 

19. Ashtekar, A., Hansen, R.O.: A unified treatment of null and spatial infinity in 
general relativity, i. universal structure, asymptotic symmetries, and conserved 
quantities at spatial infinity. J. Math. Phys. 19, 1542 (1978) 

20. Ashtekar, A.: Asymptotic structure of the gravitational field at spatial infinity. 
In: A. Held (cd.) General Relativity and Gravitation, one hundred Years after 
the Birth of Albert Einstein, vol. 2, p. 37. Plenum Press, New York (1980) 

21. Ashtekar, A., Magnon-Ashtekar, A.: On conserved quantities in general rela- 
tivity J. Math. Phys. 20, 793 (1979) 

22. Hayward, S.A.: Spatial and null infinity via advanced and retarded 
conformal factors. Phys. Rev. D 68(10), 104,015 (2003). DOI 
10.1103/PhysRevD.68.104015 

23. Frauendiener, J.: Conformal infinity. Liv. Rev. Relat. 7, 1 (2004). URL (cited 
on 8 May 2009): http://www.livingreviews.org/lrr-2004-l 

24. Arnowitt, R., Dcscr, S., Misner, C.W.: The dynamics of general relativity. In: 
L. Wittcn (cd.) Gravitation: an introduction to current research, p. 227. John 
Wiley & Witten (1962) 

25. Henneaux, M., Teitelboim, C: Quantization of gauge systems. Princeton, USA: 
Univ. Pr. (1992) 520 p (1992) 

26. Abbott, L.F., Deser, S.: Stability of Gravity with a Cosmological Constant. 
Nucl. Phys. B195, 76 (1982). DOI 10.1016/0550-3213(82)90049-9 

27. Deser, S., Tekin, B.: Energy in generic higher curvature gravity theories. Phys. 
Rev. D 67(8), 084,009 (2003). DOI 10.1103/PhysRevD.67.084009 

28. Misner, C.W.: Gravitational field energy and gOO. Phys. Rev. 130(4), 1590- 
1594 (1963). DOI 10.1103/PhysRev.l30.1590 



38 



Jose Luis Jaramillo and Eric Gourgoulhon 



29. Beig, R.: Arnowitt-deser-misner energy and p'oo- Phys. Lett. 69A, 153 (1978) 

30. Grandclcmcnt, P., Gourgoulhon, E., Bonazzola, S.: Binary black holes in cir- 
cular orbits. 11. Numerical methods and first results. Phys. Rev. D65, 044,021 
(2002). DOI 10.1103/PhysRevD. 65.044021 

31. Schoen, R., Yau, S.T.: On the positive mass conjecture in general relativity. 
Comm. Math. Phys. 65, 45 (1979) 

32. Schoen, R., Yau, S.T.: Proof of the positive mass theorem, ii. Comm. Math. 
Phys. 79, 231 (1981) 

33. Witten, E.: A new proof of the positive energy theorem. Comm. Math. Phys. 
80, 381 (1981) 

34. Deser, S., Teitelboim, C: Supergravity has positive energy. Phys. Rev. Lett. 
39(5), 249-252 (1977). DOI 10.1103/PhysRevLett. 39.249 

35. York Jr., J.W.: Kinematics and dynamics of general relativity. In: L.L. Smarr 
(ed.) Sources of Gravitational Radiation, pp. 83-126. Cambridge University 
Press, Cambridge, UK (1979) 

36. Bondi, H., van der Burg, M.G.J., Metzner, A.W.K.: Gravitational waves in 
general relativity Vll. Waves from axi-symmetric isolated systems. Proc. Roy. 
Soc. Lond. A 269, 21 (1962) 

37. Sachs, R.K.: Gravitational waves in general relativity VIII. Waves in asymp- 
totically fiat space-time. Proc. Roy. Soc. Lond. A 270, 103 (1962) 

38. Penrose, R.: Asymptotic properties of fields and space-times. Phys. Rev. Lett. 
10, 66 (1963) 

39. Wald, R.M., Zoupas, A.: A General Definition of "Conserved Quantities" in 
General Relativity and Other Theories of Gravity. Phys. Rev. D61, 084,027 
(2000). DOI 10. 1103/PhysRevD. 61.084027 

40. Geroch, R., Winicour, J.: Linkages in general relativity. J. Math. Phys. 22 
(1981) 

41. Winicour, J., Tamburino, L.: Lorentz-covariant gravitational energy- 
momentum linkages. Phys. Rev. Lett. 15(15), 601-605 (1965). DOI 
10.1103/PhysRevLett.l5.601 

42. Schoen, R., Yau, S.T.: Proof that the bondi mass is positive. Phys. Rev. Lett. 
48(6), 369-371 (1982). DOI 10.1103/PhysRevLett.48.369 

43. Horowitz, G.T., Perry, M.J.: Gravitational energy cannot become negative. 
Phys. Rev. Lett. 48(6), 371-374 (1982). DOI 10.1103/PhysRevLett.48.371 

44. Misner, C.W., Sharp, D.H.: Relativistic equations for adiabatic, spherically 
symmetric gravitational collapse. Phys. Rev. 136(2B), B571-B576 (1964). 
DOI 10.1103/PhysRev.l36.B571 

45. Brown, J.D., York, J.W.: Quasilocal energy and conserved charges derived 
from the gravitational action. Phys. Rev. D 47(4), 1407-1419 (1993). DOI 
10.1103/PhysRevD.47.1407 

46. Kijowski, J.: A simple derivation of canonical structure and quasi-local hamil- 
tonians in general relativity. Gen. Relativ. Gravit. 29, 1997 (307-343) 

47. Liu, C.C.M., Yau, S.T.: Positivity of quasilocal mass. Phys. Rev. Lett. 90(23), 
231,102 (2003). DOI 10. 1103/PhysRevLett.90. 231102 

48. Epp, R.J.: Angular momentum and an invariant quasilocal energy in 
general relativity. Phys. Rev. D 62(12), 124,018 (2000). DOI 
10.1103/PhysRevD.62.124018 

49. Wang, M.T., Yau, S.T.: Quasilocal mass in general relativity. Phys. Rev. Lett. 
102, 021,101 (2009). DOI 10.1103/PhysRevLett. 102.021101 



Mass and Angular Momentum in General Relativity 



39 



50. Murchadha, N.O., Tung, R.S., Xie, N.: Quasilocal energy in general relativity. 
arXiv : 0905 . 0647 [gr-qc] 

51. Hawking, S.: Gravitational radiation in an expanding universe. J. Math. Phys. 
9, 598-604 (1968) 

52. Huisken, G., Ilmanen, T.: The inverse mean curvature flow and the riemannian 
Penrose inequality. J. Difl^. Geom. 59, 353 (2001) 

53. Gcroch, R.: Energy extraction. Ann. N.Y. Acad. Sci. 224, 108-117 (1973) 

54. Hayward, S.A.: Quasilocal gravitational energy. Phys. Rev. D 49(2), 831-839 
(1994). DOl 10.1103/PliysRcvD.49.831 

55. Baxtnik, R.: New definition of quasilocal mass. Phys. Rev. Lett. 62(20), 2346- 
2348 (1989). DOI 10.1103/PhysRevLett.62.2346 

56. Bray, H.L., Chrusciel, P.T.: The Penrose inequality. In: P.T. Chrusciel, 
H. Friedrich (cds.) The Einstein Equations and the Largo Scale behavior of 
Gravitational Fields, p. 39. Bikhauscr, Basel (2004) 

57. Yang, J., Ma, Y.: Quasilocal energy in loop quantum gravity. Phys. Rev. D 
80(8), 084,027 (2009). DOI 10.1103/PhysRevD.80.084027 

58. Dreyer, O., Krishnan, B., Shoemaker, D., Schnetter, E.: Introduction to iso- 
lated horizons in numerical relativity. Phys. Rev. D67, 024,018 (2003) 

59. Cook, G.B., Whiting, B.F.: Approximate killing vectors on s^. Phys. Rev. 
D76, 041,501 (2007) 

60. Hayward, S.: General laws of black-hole dynamics. Phys. Rev. D 49, 6467 
(1994) 

61. Ashtekar, A., Krishnan, B.: Dynamical horizons: Energy, angular momentum, 
fluxes and balance laws. Phys. Rev. Lett. 89, 261,101 (2002) 

62. Ashtekar, A., Krishnan, B.: Isolated and dynamical horizons and their ap- 
plications. Liv. Rev. Relat. 7, 10 (2004). URL (cited on 8 May 2009): 
http: / /www.livingreviews.org/lrr- 2004^10 

63. Hayward, S.A.: Angular momentum conservation for dynamical black holes. 
Phys. Rev. D74, 104,013 (2006) 

64. Korzynski, M.: Quasi-local angular momentum of non-symmetric isolated and 
dynamical horizons from the conformal decomposition of the metric. Glass. 
Quant. Grav. 24, 5935-5944 (2007) 

65. Szabados, L.B.: On a class of 2-surface observables in general relativity. Glass. 
Quant. Grav. 23, 2291 2302 (2006) 

66. Campanelli, M., Lousto, CO., Zlochower, Y., Krishnan, B., Merritt, D.: Spin 
flips and precession in black-hole-binary mergers. Phys. Rev. D75, 064,030 
(2007) 

67. Ashtekar, A., Fairhurst, S., Krishnan, B.: Isolated horizons: Hamiltonian 
evolution and the first law. Phys. Rev. D62, 104,025 (2000). DOI 
10.1103/PhysRevD. 62. 104025 

68. Ashtekar, A., Beetle, C., Lewandowski, ,J.: Mccharucs of Rotating Isolated 
Horizons. Phys. Rev. D64, 044,016 (2001). DOI 10.1103/PhysRevD.64.044016 

69. Gourgoulhon, E., Jaramillo, J.L.: A 3-1-1 perspective on null hypersurfaces and 
isolated horizons. Phys. Rcpt. 423, 159 (2006) 

70. Crnkovic, C., Wittcn, E.: Covariant description of canonical formalism in ge- 
ometrical tlioorics. In: S. Hawking, W. Israel (eds.) Three Hundred Years of 
Gravitation. Cambridge University Press, Cambridge (1987) 

71. Christodoulou, D.: Reversible and irreversible transformations in black- 
hole physics. Phys. Rev. Lett. 25(22), 1596-1597 (1970). DOI 
10.1103/PhysRevLett.25.1596 



40 



Jose Luis Jaramillo and Eric Gourgoulhon 



72. Malec, E., O'Murchadha, N.: Trapped surfaces and the penrose inequality in 
spherically symmetric geometries. Phys. Rev. D 49(12), 6931-6934 (1994). 
DOT 10.1103/PhysRevD.49.6931 

73. Bray, H.: Proof of the ricmaniiiaii perirosc conjecture using the positive mass 
theorem. J. Diff. Geom. 59, 177 (2001) 

74. Gibbons, G.: The isoperimetric and bogomolny inequalities for black holes. 
In: T. Willmore, N. Hitchin (eds.) Global Ricmannian Geometry, pp. 194-202. 
Ellis Horwood; Halsted Press, Chichester, U.K.; New York, U.S.A. (1984) 

75. Ben-Dov, I.: The penrose inequality and apparent horizons. Phys. Rev. D70, 
124,031 (2004) 

76. Bardeen, J.M., Carter, B., Hawking, S.W.: The four laws of black hole me- 
chanics. Commun. Math. Phys. 31, 161 (1973) 

77. Hawking, S.W.; Particle Creation by Black Holes. Commun. Math. Phys. 43, 
199-220 (1975). DOT 10.1007/BF02345020 

78. Hawking, S.W., Ellis, G.F.R.: The large scale structure of space-time. Cam- 
bridge University Press (1973) 

79. Hawking, S.W.: Gravitational radiation from colliding black holes. Phys. Rev. 
Lett. 26(21), 1344-1346 (1971). DOI 10.1103/PhysRevLett. 26.1344 

80. Hawking, S.W.: Black holes in general relativity. Commun. Math. Phys. 25, 
152-166 (1972) 

81. Israel, W.: Third law of black-hole dynamics: A formulation and proof. Phys. 
Rev. Lett. 57(4), 397-399 (1986). DOI 10.1103/PhysRevLett.57.397 

82. EYiedman, J.L., Uryu, K., Shibata, M.: Thermodynamics of binary black 
holes and neutron stars. Phys. Rev. D65, 064,035 (2002). DOI 
10.1103/PhysRevD. 65.064035 

83. Booth, I., Fairhurst, S.: The first law for slowly evolving horizons. Phys. Rev. 
Lett. 92, 011,102 (2004) 

84. Gourgoulhon, E., Jaramillo, J.L.: Area evolution, bulk viscosity and entropy 
principles for dynamical horizons. Phys. Rev. D74, 087,502 (2006) 

85. Booth, I., Fairhurst, S.: Isolated, slowly evolving, and dynamical trapping hori- 
zons: geometry and mechanics from surface deformations. Phys. Rev. D75, 
084,019 (2007) 

86. Dain, S.: Angular-momentum-mass inequality for axisymmetric black holes. 
Phys. Rev. Lett. 96, 101,101 (2006) 

87. Chrusciel, P.T.: Mass and angular-momentum inequalities for axi-symmetric 
initial data sets I. Positivity of mass. Annals Phys. 323, 2566-2590 (2008). 
DOI 10.1016/j.aop.2007.12.010 

88. Chrusciel, P.T., Li, Y., Weinstein, G.: Mass and angular-momentum inequali- 
ties for axi-symmctric initial data sets. II. Angular-momentum. Annals Phys. 
323, 2591-2613 (2008). DOI 10. 1016/j.aop. 2007.12.011 

89. Dain, S., Ortiz, O.: Numerical evidences for the angular momentum-mass in- 
equality for multiple axially symmetric black holes. arXiv : 0905 . 0708 [gr-qc] 

90. Ansorg, M., Petroff, D.: Negative komar mass of single objects in regular, 
asymptotically flat spacetimes. Class. Quantum Grav. 23, LSI (2006) 

91. Ansorg, M., Pfister, H.: A universal constraint between charge and rotation 
rate for degenerate black holes surrounded by matter. Class. Quant. Grav. 25, 
035,009 (2008). DOI 10.1088/0264-9381/25/3/035009 

92. Hennig, J., Ansorg, M., Cederbaum, C: A universal inequality between the 
angular momentum and horizon area for axisymmetric and stationary black 
holes with surrounding matter. Class. Quantum Grav. 25, 162,002 (2008) 



Mass and Angular Momentum in General Relativity 



41 



93. Hennig, J., Cederbaum, C, Ansorg, M.: A universal inequality for axisymmet- 
ric and stationary black holes with surrounding matter in the einstein-maxwell 
theory. arXiv : 0812. 2811 [gr-qc] 

94. Haag, R.: Local Quantum Physics, Fields, Particles, Algebras. Texts and 
Monographs in Physics. Springer, Berlin, Germany; New York, U.S.A. (1992) 



